First-Principles Theory of Multipolar Order in Neptunium Dioxide 
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We provide a first-principle, materials-specific theory of multipolar order and superexchange in 
Np02 by means of a non-collinear local-density approximation +U (LDA+?7) method. Our calcu- 
lations offer a precise microscopic description of the triple-g-antiferro ordered phase in the absence 
of any dipolar moment. We find that, while the most common non-dipolar degrees of freedom (e.g., 
electric quadrupoles and magnetic octupoles) are active in the ordered phase, both the usually ne- 
glected higher-order multipoles (electric hexadecapoles and magnetic triakontadipoles) have at least 
an equally significant effect. 
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One of the milestone achievements in the field of 
magnetism was the discovery that the relatively sim- 
ple Heisenberg-Dirac exchange Hamiltonian successfully 
describes the intersite interactions for most 3<i-electron- 
based compounds [1]. In contrast, superexchange in- 
teraction involving terms other than magnetic dipoles 
can become effective when the orbital degrees of freedom 
are unquenched and give rise to exceptional physical be- 
haviour, such as unconventional ordered phases governed 
by multipolar order parameters [2|, |3[. Unfortunately, 
in most cases recognizing the true driving force behind 
these anomalous phase transitions is a complicated task, 
due to the large number of independent degrees of free- 
dom to be considered. For example, more than half a 
century of extensive experimental and theoretical stud- 
ies I4j- ll8| was required to reach a consensus on the real 
nature [19| of the unusual ordered state displayed by nep- 
tunium dioxide below 25 K. Np02 is now widely recog- 
nized as the archetypal compound for magnetic multipole 
ordering [2[. Part of the problem is due to the fact that, 
although crystal-field theories (and their extensions) pro- 
vided precious insight in the local symmetry of the 5 f 
levels and the possible order parameters [3|, |4J, [lg, LL7|, LL9| , 
these do not include any material-specific aspects, such 
as hybridization with the O 2p electrons, and cannot ex- 
plain why, of all materials, Np02 is special; also, mod- 
els involving a cumbersome large number of adjustable 
parameters would have to be used if all active degrees 
of freedom were to be considered; even the most recent 
theoretical advances [2|, LLj], LL9| have to rely on drastic 
assumptions [20] to simplify the calculations, notwith- 
standing the wealth of available experimental data. 

In this Letter we provide a first-principles theory of 
multipolar order in Np02- We show that non-collinear 
local density approximation +U (LDA+Z7) calculations 
can successfully overcome the above-mentioned difficul- 
ties and provide new insight in the anomalous symmetry- 
broken phase. Somewhat surprisingly, only a few first- 
principles studies [2ll423J were ever performed on Np02, 



none of these concerning multipolar ordering because 
the necessary computational framework was not devel- 
oped. Conversely, we show here that our state-of-the- 
art full-potential density-functional theory based calcu- 
lations, using the LDA+Z7 method [24] to treat the strong 
Np /-electron correlations, provide a realistic material- 
specific description of the electronic structure in the or- 
dered phase, including the role and influence of each in- 
dividual multipole. In particular, two specific properties 
of the LDA+Z7 framework make it suitable to study com- 
plex multipolar-ordered phases: the spin-orbital depen- 
dence of the local potential, which is essential since ap- 
pearing multipole order parameters involve multiple spin 
and orbital degrees of freedom, and the ability to take 
into account non-collinearity of local order parameters, 
because the multipole moments on each Np site can be 
expressed through the local potential only. 

In the one-electron scheme, the strong spin-orbit cou- 
pling splits the 5/ orbitals of a Np 4+ (5/ 3 ) ion in a j—b/2 
sextet and a j—7/2 octet. Np02 crystallizes in the cu- 
bic CaF2 (Fm3m) structure, which consists of interpen- 
etrating Np and oxygen FCC lattices. Under a cubic 
crystal field the lower j=5/2 multiplet splits into a T$ 
quartet and a T^ doublet. In the transition to the mul- 
tipolar ordered state the atomic distances in the crystal 
structure are unchanged, yet the Np site-local cubic Oh 
symmetry is broken. The (111) orientation of the 3q- 
antiferro electric quadrupoles, detected in resonant x-ray 
scattering (RXS) experiments [5j, implies a lowering of 
the structural symmetry to that of a 4-Np02 sublattices 
SC unit cell, where the Np site has local Dsd symme- 
try [25[. As a result, the j=5/2 states are split in two 
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doublets, Tg ; and T K 6 , and two singlets, T4 and T^; the 
latter two are actually degenerate if time-reversal sym- 
metry is not broken [26[. Our selfconsistent calculations 
were performed using the full-potential linearized aug- 
mented plane wave (FLAPW) method, as implemented 
in the TSPACE and KANSAI program packages, com- 
bined with the LDA+U approach [27] to treat properly 
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FIG. 1: (Color online) Computed charge density cross- 
sections of Np02- a,b, Non- multipolar ordered state with 
cubic FCC symmetry in the [100], respectively, [110] plane. 
c,d, Magnetic multipolar state showing the (111) 3q antiferro 
electric long-range order in the [100], respectively, [110] plane. 
Both solutions were obtained with U=4 eV and J=0 eV. 
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FIG. 2: (Color online) a, Orbital-projected density of states 
(DOS) of a multipolar state computed with [7=4 eV, J=0 
eV. Plotted is the Np 5/, j-=5/2, DOS projected on compo- 
nents of the irreducible Dsd basis, and the Np j=7/2 and O 
p character DOS. b, Orbital-projected DOS of the multipo- 
lar phase obtained with U—A eV and J=0.2 eV. c, Calcu- 
lated LDA+[7 one-electron orbital occupation numbers and 
off-diagonal (T Q , Tq ^) matrix elements as a function of ex- 
change J at [7=4 eV. d, Calculated expectation values of the 
multipolar (quadrupolar, octupolar, hexadecapolar, and tri- 
akontadipolar) order parameters in Np02 as a function of J. 



the on-site Coulomb correlations of the Np / electrons. 
The relativistic spin-orbit interaction as well as the lo- 
cal +U potential were included via a second variational 
step. For the LDA+/7 double counting term the fully 
localized limit was adopted. Full-potential LDA+/7 was 
shown to describe an anisotropic /-state [27]] , which is es- 
pecially important for the multipole ordered state. The 
above-mentioned unit cell and the local Dsd symmetry 
at the Np sites, together with broken time-reversal sym- 
metry to accommodate magnetic solutions, have been 
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used. Our Z^d-irreducible basis states Tq (r 4 , T 6 
r$) are constructed from the T? (Tg) orbitals in Oh 
symmetry. For the on-site Coulomb U and exchange J 
we use values of 3-5 eV, respectively, 0-0.5 eV, which 
are in the commonly accepted range for light actinides. 
The two-electron integrals of the Coulomb interaction 
of the / electrons were expressed in terms of effective 
Slater integrals F K (0<ft<6), which are associated with 
F =U, and, F 4 =(41/297)F 2 , F 6 =(175/11583)F 2 , and 
J=(286F 2 +195F 4 +250F 6 )/6435 @. 

Fig. [T] shows calculated cross-sections of the charge 
density of Np0 2 in the high-symmetry [100] and [110] 
planes. Panels a and b depict cross-sections for a non- 
multipolar ordered state, obtained for the FCC unit cell. 
Panels c and d show cross-sections for a magnetic multi- 
polar ordered phase, computed for the 4-sublattices SC 
unit cell, with the same values of the Coulomb U and ex- 
change J parameters. A comparison of the charge densi- 



ties in panels a, c and b, d reveals the symmetry reduction 
of the charge distribution in the multipolar ordered state, 
which is found to have a significantly deeper total energy 
than that of the non-magnetic single FCC unit cell (in 
which multipolar order cannot occur). The 3qr-antiferro 
electric long-range order can be recognized clearly from 
panel d, where the charge lobes are tilted and point along 
(111) directions. 

Fig. |2H ? b show the partial density of states (DOS) of 
the multipolar state computed with U=A eV, J=0 and 
U=4 eV, J=0.2 eV, respectively. The Np 5/, j=5/2, 
DOS is projected on components of the irreducible Dsd 
basis. The computed DOS reveals that Np0 2 in the 
multipolar ordered phase is an insulator with a gap of 
2 eV (£7=4 eV, J=0 eV). In sharp contrast, cubic non- 
magnetic Np0 2 is computed to be metallic. The orbital 
projected DOS panels a and b evidence a significant hy- 
bridization between oxygen 2p orbitals and Np 5/, j=5/2 
orbitals. The tiny contribution of the Np j=7/2 states 
to the valence states provides a ground to neglect them 
as was assumed in mean- field models [17[. Fig. [2fc shows 
the calculated diagonal (occupation number) and non- 
diagonal (r^ , Tq ) density matrix element as a function 
of J, with U fixed at 4 eV. When J=0, the occupied Np 

5/ orbitals have strong T4 and T 6 ' character whereas 

(2) 
the unoccupied orbitals have T$ and T 6 character. As 



soon as a non-zero value of J is taken into account, the 
orbital character changes due to a mixing of the Ig ^ and 
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Tq states. The occupation of the Tq orbitals starts to 
increase significantly, at the expense of the Tq ) orbitals, 
and a relatively large off-diagonal term between the two 
doublets is present, which saturates for J>0.15 eV. In 
order to obtain microscopic insight into the mechanisms 
which drive the phase transition, we calculated the expec- 
tation values for the different multipoles having T$ sym- 
metry and ranks from 2 to 5 (electric quadrupoles, mag- 
netic octupoles, electric hexadecapoles, and magnetic tri- 
akontadipoles) as a function of the exchange J (Fig.|2]i), 
using their tensor-operator expressions OlL> [29|, |30j: 
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where \kb) are the Bloch states, r labels the atomic site, 
£ the angular momentum quantum number, and 7 (7') 
indexes combined spin and orbital degrees. The latter 
can be expressed in an £, m, s basis, or via a unitary 
transformation, in the irreducible basis. nlL, is the den- 
sity matrix in the muffin- tin sphere. Each multipolar 
moment has been normalized to the value which would 
be obtained if only the three lowest orbitals in the one- 
electron description, i.e., the T4 singlet and T Q doublet, 
were filled. Our ab initio study demonstrates that, while 
the occupation of the j=7/2 octet can be safely neglected, 
all the orbitals which arise the j=5/2 sextet largely con- 
tribute to the ordered phase, with the effect to boost the 
expectation values of the generally discarded higher-rank 
multipole moments. 

We consider first the magnetic multipoles, whose the 
primary order parameters belongs. For J=0, their 
normalized expectation values coincide, an indication 
that triakontadipoles are inactive (and undistinguishable 
from octupoles). However, for non-zero J values the 
(Tg ,T 6 ^) off-diagonal term increases. This mixing "un- 
locks" the triakontadipolar degrees of freedom; whilst the 
expectation value of the octupoles is significantly reduced 
(changing its sign around J=0.04 eV and saturating at 
only about half of the original value for J>0.15 eV), the 
triakontadipolar moment is strongly enhanced. This find- 
ing is compatible with previous model calculations [19| 
yielding a very low octupolar polarizability for the or- 
dered ground state and unveils the microscopic mecha- 
nism which leads the T$ triakontadipole to become the 
primary order parameter. 

The effect observed on the electric multipoles as sec- 
ondary order parameters is even more impressive. The 
influence of the Tq -Tq coupling is so large for the 
quadrupole that the normalized expectation values of the 
two electric multipole moments are different in magni- 
tude and sign already when J=0 eV. Increasing the J 
value leads to a reduction of the hexadecapole and to a 
change in sign of the quadrupole, but the former remains 




FIG. 3: (Color online) Computed charge and spin distribu- 
tions of the Np 5/ electrons in NpC^. a, Non- multipolar 
ordered state, computed with U—A and J=0 eV. b, Magnetic 
multipolar ordered state, for [7=4 and J=0 eV. c, Magnetic 
multipole state, for U—A and J=0.04 eV. b, Magnetic multi- 
pole state, for [7=4 and J=0.2 eV. The charge distribution is 
depicted by the three-dimensional hypersurface and the direc- 
tion of the current flow is shown by the thin lines. The colors 
on the three-dimensional charge distribution depict the local 
spin magnetization projected on the z quantization axis. Red, 
blue colors label, respectively, a net up or down spin. Green 
color labels a vanishing net spin magnetization. 



active and well-distinguished from the latter in the whole 
investigated J range. On these grounds, it might seem 
quite striking that the role of Y$ hexadecapolar degrees 
of freedom in the ordered phase of Np02 has never been 
considered so far, but it must be reminded that hexade- 
capolar order could not be detected by a direct RXS ex- 
periment (in contrast to quadrupoles) due to the very 
small cross-section of E2-E2 transitions [5[ . 

Fig. [3] shows calculated charge and magnetization dis- 
tributions of the 5/ electrons on a Np ion, for U=4 eV 
and several J values. Fig. [3^ depicts the charge and spin 
distribution of the non-ordered state. The cubic T$ sym- 
metry of this distribution can be seen from the charge 
lobes pointing along the cubic axes of the CaF2 struc- 
ture. The green color exemplifies that the distribution is 
non-magnetic everywhere. Panel [3)3 shows the distribu- 
tion of the magnetic multipolar ordered state (for J=0 
eV). The local symmetry axis of the charge distribution 
has become the (111) axis. Magnetic multipolar order is 
revealed by the pattern of the self-consistent spin distri- 
bution: a net up (down) magnetization along the z axis 
is depicted by the red (blue) color. The dominantly red 



color occurring on the charge hypersurface around the 
(111) direction indicates a net up-magnetization along 
the (111) axis. Fig. [3h shows the multipolar distribu- 
tion computed for J=0.04 eV; an inspection of Fig. [2] 
shows that for this J value the distribution contains no 
quadrupolar nor octupolar contribution, only a hexade- 
capolar and triakontadipolar contribution. The fact that 
this multipole looks similar to the one in panel [3]o allows 
one to infer that when the full j=5/2 multiplet is consid- 
ered in selfconsistent calculations, in contrast with pre- 
vious results obtained by simplified models, these high- 
rank contributions are always important and at least as 
significant as lower-rank ones; in particular, the shape 
of the charge distribution in the small-exchange region 
is mainly due the electric hexadecapoles, whilst mag- 
netic triakontadipoles dominate over the octupoles for 
all reasonable choices of J. Fig. [SJi presents the ob- 
tained magnetic multipolar distribution for J=0.2 eV. 
This charge distribution maintains the symmetry along 
the (111) axis, but is more spherical. This can be under- 
stood as the result of a competition between the multipo- 
lar interaction and the crystal-field potential; the former 

(2) 

tends to increase the occupancy of the T K 6 ' doublet to 
the expense of states preferred by the latter, so that the 
cubic symmetry of the lattice is less evident (equal occu- 
pation of all states belonging to the whole j=5/2 sextet 
would automatically lead to spherical symmetry). The 
magnetic multipoles computed for larger J values bear 
no significant difference from the one for J=0.2 eV. 

Our state-of-the-art calculations offer new insight in 
the exceptional symmetry-broken phase of NpC>2. The 
local symmetry breaking resulting from the ordering of 
both electric and magnetic multipoles, with the excep- 
tion of dipoles, significantly lowers the total energy with 
respect to the non-magnetic state. Multipolar superex- 
change contributions drive Np02 to an insulating 3q- 
ordered ground state. The anomalous phase originates 
from a strong coupling of spin and orbital degrees caused 
by the relativistic spin-orbit interaction and a compe- 
tition between the crystal-field potential and the mul- 
tipolar superexchange. Whilst the most common non- 
dipolar degrees of freedom (e.g. electric quadrupoles and 
magnetic octupoles) are present and active, we find that 
both the usually neglected higher-order multipoles (elec- 
tric hexadecapoles and magnetic triakontadipoles) play 
at least an equally important role. Our calculations em- 
phasize the important influence of the higher- lying states 
in the multipole formation, and reveal that a strong mix- 
ing between Tg ^ and Tg ^ states, due to an increased ex- 
change J, causes a steep growth of the triakontadipole at 
the expense of the octupole. Putting our first-principles 
results in context with the available experimental data, 
this allows to infer that the primary order parameter is 
the T 5 triakontadipole. 
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